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Abstract 

We prove the stability of a spacetime with a naked singularity under scalar 
field perturbations, where the perturbations are regular at the singularity. 
This spacetime, found by Janis, Newman and Winicour, and independently 
by Wyman, is sourced by a massless scalar field and also arises as a certain 
limit of a class of charged dilatonic solutions in string theory. This stability 
result opens up specific questions for investigation related to the cosmic cen- 
sorship conjecture and the mechanism by which it is implemented in nature. 
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I Introduction 



There is a vast body of work on the issue of classical stability of black holes 
under perturbations. As is well-known, the region exterior to the horizon of 
the Schwarzschild black hole is linearly stable under perturbations by various 
fields [lj. Specifically, there are stability results for perturbations (for e.g., 
scalar field or gravitational perturbations) that are initially of compact sup- 
port and obey Dirichlet or Neumann boundary conditions at the horizon and 
infinity. These stability results have also been shown in Kruskal coordinates, 
which are valid across the horizon. However, since the perturbations live in 
the spacetime domain exterior to the black hole horizon, they cannot access 
the black hole curvature singularity, and it is not the boundary at which 
boundary conditions are enforced. □ 

However, one could now ask similar stability questions for situations in 
which the perturbations live in a spacetime domain with a curvature sin- 
gularity, or when the curvature singularity is a natural boundary at which 
boundary conditions need to be enforced. An example of this is the question 
of stability of a spacetime containing a naked singularity under various types 
of perturbations. Such a spacetime would be geodesically incomplete, as can 
be probed by point particle probes which move along geodesies. Suppose, 
instead of point particle probes, if fields were used as probes (i.e., the space- 
time was perturbed by a field), then would the time evolution of the field 
be well-defined in a spacetime with a naked singularity? This question has 
been addressed by Wald [2] and later by many authors (see, for example, 
[2], [I], [5], [6] and references therein). In [7], the authors explicitly show an 
example where the time evolution of various fields is well-defined in such a 
spacetime. If the time evolution of fields is well-defined, then the next ques- 
tion is one of the stability of the spacetime. If the spacetime is perturbed 
by a field configuration which is initially of compact support, for example, 
then how does this perturbation evolve? The stability of the negative mass 
Schwarzschild spacetime (whose curvature singularity now becomes naked, 
as there is no horizon) under gravitational perturbations was analyzed in jS] 
- with boundary conditions demanding regularity of the perturbation at 
the singularity, this spacetime was shown to be unstable [9]. The Reissner- 



3 Other classes of perturbations, such as quasinormal modes have ingoing boundary 
conditions at the horizon, which are motivated by the presence of the curvature singularity 
in the black hole interior. However, they are not normalizable and are thus different from 
the perturbations considered in the usual classical stability analysis. 
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Nordstrom metric with charge larger than mass (super-extremal) and the 
Kerr metric with angular momentum larger than mass are spacetimes with 
naked singularities. They were shown to be classically unstable [ID]. It was 
shown by Christodoulou [TT] that a naked singularity formed during self- 
similar collapse of a spherically symmetric massless scalar field was unstable 
(a similar instability result can be found in [12])- These instability results 
seem consistent with cosmic censorship and the view that naked singularities 
are not likely to form from generic initial conditions in collapse situations. 
One example of a naked singularity stable under gravitational perturbations 
is the self-similar Vaidya spacetime which is sourced by dust. However this 
is not a serious challenge to cosmic censorship since the matter sourcing 
the spacetime is null dust, which forms singularities under collapse in flat 
spacetime [13]. 

In this paper, we probe a naked singularity spacetime which is a solution 
to the Einstein equation with spherical symmetry and sourced by a mass- 
less scalar field. The spacetime, which we call the Janis-Newman-Winicour 
(JNW) metric, was found independently by Janis, Newman and Winicour 
[13] and Wyman [13] who used different coordinates from each other — the 
two solutions were later shown to be identical, and containing a naked singu- 
larity [IE] , [TT] . This spacetime can also be obtained from a class of charged 
dilaton black holes found by Gibbons and Maeda which are described by 
mass, charge and a coupling strength of the dilaton to the Maxwell field 
[T5] (the solutions and some features were discussed independently later by 
Garfinkle, Horowitz and Strominger [E]). The JNW spacetime is obtained 
by setting the charge to zero in this class of spacetimes (resulting in a naked 
singularity), and the dilaton plays the role of the scalar field sourcing the 
JNW spacetime. More precisely, the JNW spacetime is obtained from the 
charged dilaton metrics in [19J by setting r + = 0, with r_ ^ in equations 
(20-23) of that paper. For other values of r + and r_ which correspond to 
charged dilatonic black holes, evolution of scalar fields in such spacetimes are 
discussed in (3]. Also, we note that in [5], the scalar field equation in the 
background of the JNW solution (in the coordinates of Wyman) is stated 
as an example that is well-posed (when the function space is an appropriate 
Sobolev space), making the spacetime 'wave-regular'. However analysis of 
its stability has not been carried out in that paper. 

When the magnitude of the scalar field is zero, the JNW metric reduces 
to the Schwarzschild black hole metric, and for any nonzero magnitude of the 
scalar field, the event horizon of the black hole now degenerates into a naked 
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singularity. We probe this spacetime with a test scalar field and find that 
the JNW spacetime is stable with respect to scalar field perturbations - 
this result is irrespective of the mass of the probe field. This is in contrast to 
studies of perturbations of various naked singularities, like the super-extremal 
Reissner-Nordstrom metric and Kerr metric with large angular momentum, 
which are unstable (to gravitational perturbations). Of course, it remains 
to be seen whether the JNW spacetime continues to be stable with respect 
to gravitational field perturbations — this is work in progress. If this were 
true, it would raise the question of whether the JNW spacetime can form in 
a realistic situation of collapse of a scalar field for generic initial conditions, 
and whether this stability result poses a challenge for cosmic censorship. 

In the next section, we introduce the JNW metric, and rigorously prove its 
stability under scalar field perturbations. We do this by first putting appro- 
priate boundary conditions on the perturbation at the naked singularity, and 
then showing that a perturbation that is initially bounded remains bounded 
for all time. We conclude with a section discussing the result and the in- 
teresting questions that arise from it, particularly relating to whether the 
cosmic censorship conjecture needs modification, and how cosmic censorship 
is implemented in nature. 



II Scalar field in naked singularity 

The Janis-Newman-Winicour (JNW) metric, which is a solution to the Ein- 
stein equation with spherical symmetry in the presence of a massless scalar 
field [2], [H] is given by 



ds 2 = -(1 - bjr) v dt 2 + _\ /r y dr2 + r2 (! - b/rf- v dn 2 ] (II.l) 

and dVt 2 is the standard metric on a unit two-sphere. 
The scalar field which is a source for this spacetime is 

$ = — ^=ln(l-6/r). (II.2) 
by An 



The parameter b appearing in the JNW metric is related to the ADM mass 
of the spacetime M by b = 2-y/ M 2 + q 2 . v = and when u — 1 (q — 0), we 
recover the Schwarzschild metric as there is no scalar field. As we increase q 
(the magnitude of the scalar field), v decreases from 1 to 0. An inspection of 
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the curvature invariants reveals that there is a singularity at r = b when < 
v < 1 for the Ricci and Weyl scalars as well as the Kretschmann invariant, 
and we will assume this range for v throughout the paper. This singularity is 
globally naked and the weak energy condition is satisfied for this spacetime 
[16] , [17] . Gravitational lensing near this singularity was studied in [20] , and 
later in [21] . Further properties of point particles in the JNW spacetime and 
their astrophysical consequences were studied in [22J and [23]. It is evident 
from these studies that there are pathologies both in the behaviour of massive 
particles and light rays in this spacetime, and clear differences between their 
behaviour near the naked singularity and the behaviour of particles near the 
Schwarzschild horizon. 

The question we want to investigate is that of the stability of this naked 
singularity spacetime under perturbations. The simplest perturbation we 
can consider is by another test scalar field of mass m, whose magnitude is 
everywhere smaller than the source scalar field (so that we can use linearized 
perturbation theory). Let us denote this test scalar field \. Assuming that 
the only interaction of this scalar field is through minimal coupling with the 
background gravitational field, its action is 



S = \/—ad x 



(II.3) 



The equation of motion of the test scalar field \ arising from this action 
is given by the massive Klein-Gordon equation, 

1 ^[^g^d KX }=m 2 x- (II.4) 



where g denotes the determinant of the metric (IH.lj) . 



Let f(r) = (1 - b/rf. 



fin 



We now choose the following ansatz for x : 

x =±Y lm (Q^)e^ (H.5) 
r 

Yi m (6,(j)) are the spherical harmonics, and u can be complex. We will 
prove in two stages that with reasonable physical boundary conditions, the 
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only solutions are those with oo real. This implies that there are no solutions 
growing exponentially in time. 

With the ansatz (111.5j) . equation fill. 41) becomes 



d 
dr 



'1 _ b ^ 
r dr 



+ 



b 1(1 + 1) 2 1 
s + n +m 2 - 



f 



b\ {1 - 2u) 
1-- i,. (II 

r / 



Define the coordinate r* by dr* = dr/(l — -). Note that i/iis is not t/ie 
wswa/ tortoise coordinate that changes the r — t part of the metric into a 
conformally flat form. This coordinate change is mainly for convenience in 
analyzing f ]II.6j) . Then ( 1II.6j) becomes 



d 2 jj 
dr 2 



+ (1- 



1(1 + 1) 



'I - 



m 



f 



1> 



ij. (11.7) 



r* = r + 61n(r — b), and so, as 6 < r < oo, — oo < r* < oo. The range of 
r* is infinite; however, for the actual tortoise coordinate of this metric, the 
location of the singularity is at a finite value of the tortoise coordinate. Due 
to the range of the new coordinate r* not being finite, ( jTT7j) is a singular 
Sturm-Liouville differential equation of the form 



Mij) = oj 2 w(r)^p, 



(11.8) 



where 



M 



w(r) 



d 2 6, 



, !(* + !) , a(l- 
H ^ h m 



6 

1 - - 

r 



2(1-0 



./ 



(II.9) 



First we would like to prove that u 2 has to be real. This follows if M, act- 
ing on the space of functions which are square integrable with respect to the 
measure w(r)dr* is self-adjoint. For this singular Sturm-Liouville problem, 
the relevant question is whether there exist a choice of boundary conditions 
for which M is self-adjoint (or has a self-adjoint extension). @ Our function 



4 In [5], it is mentioned that M acting on an appropriate Sobolev space, is essentially 
self-adjoint — however, no further details are given. 
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space is the space of square integrable functions with respect to the measure 
w(r)dr*, which is the natural choice for the Sturm-Liouville problem (III.9p . 
Therefore, we would like to do a careful analysis of self-adjointness of M on 
this function space and find the appropriate physical boundary conditions 
to impose on the perturbation. For this, we follow the standard procedure 
in literature on ordinary differential equations (see for example, the classic 
paper of Krall and Zettl [21], or the standard book [22]) and classify the 
endpoints of the interval on which the equation is defined, as regular or sin- 
gular (of limit point type or limit circle type). The endpoints in this problem 
are r* = ±00. They are not finite, hence singular. For the Sturm-Liouville 
problem defined on an interval (—00, 00) with (u> 2 ) G C (complex numbers), 
the singular endpoint —00 is of limit circle type if all solutions to (1II.9P are 
square-integrable with respect to the measure w{r)dr* on (—00, (3) for some 
G (—00, 00). This classification is independent of u 2 — if this statement is 
true for some (u 2 ) G C, then it is true for all (a; 2 ) G C (for a proof, see [26]). 
A similar definition holds for the other endpoint. If a singular endpoint is not 
of limit circle type, it is of limit point type. Note that the endpoint +00 being 
of limit point type means there must exist at least one (a; 2 ) G C for which 
some solution to ( 1H.9|) is not square-integrable with respect to the measure 



w(r)<ir* on (/3, 00). For endpoints of limit circle type, a careful choice of 
boundary conditions is needed to ensure self-adjointness of M. This analysis 
of self-adjointness is equivalent to the approach in terms of deficiency indices 
that is more commonly found in quantum mechanics books [21] • We will 
summarize the choice of boundary conditions that lead to self-adjointness, 
but let us investigate first whether the endpoints in our problem are of limit 
circle type. 

(i) Endpoint r* = —00: 



As r* — > —00 (that is, r — > b), r* ~ 61n(r — b). In this limit, (1H.7|) 
becomes 

drl b 

° = h + l1 ^- ™ 



Consider r* — >■ — 00 and < v < 1/2. We can then replace flll.lOj) by the 



7 



further approximate equation 

-^ + ?e r,/ V = 0. (11.11) 

We change variables to yfx, where x = je r *' b and rewrite the above as a 
modified Bessel equation whose solutions are the modified Bessel functions 
of order 0. As r* — > — oo, the general solution is of the usual form 

il> ~ A 1 I (2b^e r *' 2b ) + A 2 K (2b^e r *' 2b ). 

In the limit r* = — oo, what is relevant is the behaviour of the modified Bessel 
functions for small argument. In this limit, 



I {2b^e r * /2b ) ~ const.; K {2b^e r * /2b ) ~ (const.)r, 



Both these linearly independent solutions are square-integrable with respect 
to the measure w(r)dr* on (— oo,/3) for some /3 G (—00,00), because of the 
exponential fall-off of the measure function w(r) as a function of r* in the 
limit r* — > —00. Note that this statement is independent of u 2 . Thus this 
endpoint is of limit circle type when < v < 1/2. 
Consider r* — > —00 and v — 1/2. Then (III. 101) becomes 

- :rW + (^-^)e r * / V = 0. (11.12) 

arf b 

This equation is similar to (III. lip except for replacing ~ in the previous case 
by ( ^ ) G C. As can be checked, the statements made in the previous case 
are valid in this well, and the endpoint is of limit circle type. 

For 1/2 < v < 1, the equation approximating (III.lOp further is 

- iU - wV 2 * 1 ""^-^-^ = 0. (11.13) 
dri 

We now change to the (complex) variable y/x, where x = b 2(f- v) e 2 ( 1 ~' y ) r '*/ b . 
We then again get a linear combination of modified Bessel functions as a 
solution. The argument of these functions differs from the case v < 1/2, but 
similar arguments ensure that the endpoint is of limit circle type, independent 
of u 2 . 
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So we conclude that the endpoint r* = — oo is of limit circle type for 
all < v < 1 (also, the result is independent of the scalar field mass m - 
the m-dependent term is negligible near this endpoint). This means that we 
need to choose our boundary conditions carefully at this endpoint. 

(ii) Endpoint r* = +00: 



In the limit r* — > 00 (that is, r — > 00), (1II.7j) becomes 

d 2 

- — ip + m 2 tjj = u 2 ip. (H.14) 

The general solution in this limit is 

</> ~ B ie ikr * + B 2 e- ikr *- 

where k 2 = (u> 2 — m 2 ). Clearly there exist choices of (u 2 ) G C for which k 
has a nonzero imaginary part (for e.g., u 2 real and satisfying u> 2 < m 2 ). In 
this case, one of the two linearly independent solutions, e tkr * or e~ lkr * will go 
exponentially to zero as r* — > 00, and the other would grow exponentially. 
The measure w(r)dr* ~ dr* in this limit, and thus, only one of these two 
solutions is square integrable with respect to this measure. We conclude that 
the endpoint r* = +00 is of limit point type. 

Now, in order to make this Sturm-Liouville problem — i.e., the operator 
M defined in ( 111.91) — self-adjoint, we need to place appropriate boundary 



conditions which are dictated by this classification. We will closely follow 
the paper of Krall and Zettl [21] for this (including their notation). 
Some notation: we denote the endpoint a := —00, d := +00 . Then a is of 
limit circle type and d is of limit point type. Also denote 

ip(a) = lim ip{r*) = lim ip{r*). 

We can define if)(d) similarly. 

We need no boundary conditions at endpoint d since it is of limit point 
type and one of the two linearly independent solutions near d is excluded 
anyway as it is not square integrable. To specify the self-adjoint boundary 
conditions at a, we first look for solutions to the equation My = 0, i.e 



dri r 



b 1(1 + 1) 2 (1- 
-5 +^^ + m — J 



y = 0. (11.15) 
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Denoting the Wronskian of two functions / and g by W(f, g) = fg' — gf, 
we need to choose two solutions to the equation fill. 1 5j) . y\ and y 2 so that for 
the entire range of r*, 

W( yi ,y 2 ) = l. 

We note that standard results (for e.g, Abel's identity) imply that for a 
homogeneous second order ODE of the type (111.150 . the Wronskian of two 
solutions remains constant in the entire range of r*. Therefore, we can look at 
flII.15j) in the limit r* — > — oo, and identify two linearly independent solutions 



yi and y 2 with W(yi,y 2 ) = 1 in this limit. But the approximate equation in 
this limit is precisely fill. lip . Therefore, two linearly independent solutions 

(in this limit) are y 1 = A l I Q (2b^e r ^ 2b ) and y 2 = A 2 K (2b^e r ^ 2b ). 
Recall that the Wronskian of the two modified Bessel functions Io(z) and 
K (z) with respect to the variable z is (const.) jz. Since z = 2b\J^c r *l 2b , 
the Wronskian of this pair of modified Bessel functions with respect to the 
variable r* is a constant. By choosing A\ and A 2 appropriately, we can find 
yi and y 2 so that their Wronskian is unity. Thus, we have shown that there 
exist two solutions to (III. 151) . y\ and y 2 with W (1/1,2/2) = 1- 

We now return to the problem of finding self-adjoint boundary conditions 
at the endpoint a = —00. The correct boundary conditions are given by 

ciW(V>, Vi) (a) + c 2 W(ip,y 2 ) (a) 

= ci lim [ipy[ - yxif)'] + c 2 lim [ipy' 2 - y 2 if)'] = 0. (11.16) 

Any constants C\ and c 2 are allowed, provided at least one of them is 
nonzero. Our notation in the previous equation was that 

W(ij, yi )(a)= \imW(i>, yi ); W(i>,y 2 )(a)= lim y 2 ). 

As — > a + , yx — > const, and y[ — > (prime denotes derivative with respect 
to r*). As r* — > a + , y 2 — > —00 and y 2 — > const. Now we know from the 
computations done at the beginning of this section that as r* — > a + , ip is in 
general, some linear combination of the two modified Bessel functions whose 
argument depends on the value of the parameter v. For example, when 
< v < 1/2, 



Ai I (2b^c^ 2b ) + A 2 K (2b^c^ 2b )- 
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In this limit, the Bessel function K is singular, whereas Iq goes to a constant. 
Clearly, from the general form of ip and the behaviour of y\ and y% as — > 
a + , there exist some choices of c\ and C2 for which the resulting boundary 
conditions flII.16j) imply that ip is proportional to the Bessel function Iq (i.e., 



contains no term proportional to K ). Such boundary conditions would result 
in a regular if) since K Q is not regular in this limit. For example, consider 
Ci = 1 and C2 = 0. We then obtain the boundary condition 

ipy'i(a) -il>'vi(a,) = 0. 

Studying the behaviour of solutions as r* — > a + „ we see that this implies 
if}' (a) = 0. This boundary condition picks the modified Bessel function J 
(whose derivative with respect to is zero at a) and excludes K Q (whose 
derivative with respect to r* is nonzero at a). It is also easy to see that if 
C2 7^ 0, then if; will have a component that is proportional to Kq and will 
thus not be regular at a. 

We now argue that from physical considerations, we should choose the 
boundary condition resulting in a if} that is proportional to Iq (with no con- 
tribution from the K piece). We would like this scalar field to be a perturba- 
tion, which means its back-reaction can be neglected. Recall that the energy 
momentum tensor of a scalar field x °f mass m is 

Tfj, v = d.xduX ~ \gA9 a ^ a xd P X + m 2 X 2 ). (11.17) 

Consider the trace of the energy momentum tensor T = g^ v T^ v — for x 
given in terms of ip by the ansatz flll.5j) . we can verify that if ip contains a 



component proportional to Kq as r — > b, then T diverges at r = b, and thus 
this scalar field can no longer be considered a 'perturbation'. Thus the only 
physical boundary condition is that with C\ nonzero and c<i = in flII.16p . 
and without loss of generality, we can take c\ — 1. In this problem, this 
boundary condition turns out to be equivalent to the 'regularity boundary 
condition' discussed, for example, in the context of naked singularities in [7]. 

We close with a comment on the various details of our computation. To 
make our analysis easier, we studied (1II.7P which was obtained from (1II.6P 
by a coordinate change from r to r* and writing the differential equation 
in Sturm-Liouville form. There is a natural measure associated with such 
a problem, w(r), and this determines our function space (in this case, all 
square integrable functions with respect to this measure). The natural ques- 
tion that arises is, what would happen to the main features of the analysis, 
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if we had done a different coordinate change, say from r to the genuine tor- 
toise coordinate (which would change the measure). It can be seen that the 
main features of our analysis — an endpoint of limit circle type, and one 
unambiguous self-adjoint boundary condition corresponding to regularity of 
the perturbation all remain in the analysis of the resulting Sturm-Liouville 
problem. 

Our conclusion is that, with this physically motivated boundary condi- 
tion, the Sturm-Liouville problem (1II.7j) defined on the space of square in- 
tegrable functions with respect to the measure w{r)dr* is self-adjoint, and 
therefore its spectrum (given by values of u 2 ) is real. However, we will not be 
guaranteed stability unless the self adjoint operator is positive, i.e, u 2 > 0, 
so that u is real. We will now attempt to show this. 

Proof that for normalizable solutions, u is real: 

We have already shown that u 2 is real. Let us assume that u 2 < 0, so 



that oj is imaginary. Denote 



V 



b_ + ijj_ 



-u 2 > 0. Let 
J1-* 



+ m 



f 



Q 2 1 



Then we can rewrite ( 111.71) as 



d 2 j) 
dr 2 



+ Vip = 



(11.18) 



We have assumed that Q 2 > 0. So V > 0. With this choice of potential 
V, consider the Schrodinger equation 



d 2 *P 
dr 2 



+ Vi[> = Eip. 



(11.19) 



Then the equation (III. 181) corresponds to the zero eigenvalue equation 
(E = 0) in ( 111.191) . Solutions to ( III. 191) with E ^ are not connected to our 
problem, but for the moment, let us consider the eigenvalue problem for the 
Schrodinger equation ( III. 191) . Since V > 0, this implies that for normalizable 
eigenf unctions to (1II.19j) . E cannot be negative. To see this, we multiply 
(1II.19|) by ip (complex conjugate of ip) and integrate over the entire range of 
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r*. Then after integration by parts, we see that 



-ijj— 
dr* 



+00 



+ 



dip 



dr* 



dn + / V\ib\ 2 dn 



E 



2 dr*. 



(11.20) 



Now we impose the same boundary conditions on ip that led to self-adjointness 
for operator M. This ensures that the boundary term in (1II.20j) goes to zero. 
Then it is clear that the lefthand side of (111.20p is nonnegative due to V being 
nonnegative. Thus the right hand side cannot be negative, and E > 0. So 
if E = is an eigenvalue of flll.20j) . it is the ground state eigenvalue. We 
are interested in the ground state eigenfunction which we denote ipo — this 
solves (III . 1 8 [) and we would like to know if such a normalizable eigenfunction 
exists. Recall that the ground state eigenfunction for a Schrodinger problem 
has no nodes (zeroes)in the interval (—00, 00) (for a proof, see Since ip 

is continuous, we can therefore assume without loss of generality that ipQ > 
in (—00, 00) (and if it approaches zero, it only does so asymptotically). Con- 
sider the equation obeyed by ipo' 



d 2 ip 
dr 2 



+ Vipo = 0. 



(11.21) 



Then 



d 2 ipc 
dr 2 



Vip > 0, 



since we can assume ipo > 0, and V > 0. Therefore -4^- is an increasing 
function of r*. We already know the behaviour of i/)o as r» — > —00 from 
the discussion on this endpoint being of limit circle type. For the choice of 
boundary conditions made, ip is proportional to a modified Bessel function I 
and as can be checked from the properties of these Bessel functions, ^ > 

for r* very negative. Therefore, ^ remains positive. This implies ipo is an 
increasing function of r*, and furthermore, the modified Bessel function I — > 
1 as r* — t- —00. It is therefore impossible for ip to go to zero asymptotically 
as r —7- 00 as required by the boundary conditions — to do so, it must 
decrease somewhere, but its slope is never negative. Thus we conclude that 
there is no normalizable solution to the equation (1II.21j) . Our argument relied 
on the positivity of V, which in turn is due to assuming that u> 2 < 0. The 
conclusion therefore is that there are no normalizable solutions to (\ll.7h for 
u zero or pure imaginary. 
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These arguments show that with the appropriate boundary conditions, 
M is a positive self-adjoint operator acting on the Hilbert space of square 
integrable functions with respect to the measure w(r)dr* which we denote 
by H. Rewriting the Klein-Gordon equation flll.4j) as 



^ = — ( IL22 ) 

at 1 w(r) 

we start with initial data of compact support, denoted by i[>t=o and (d t ip)t=o, 
and investigate if the solutions to (1II.22j) that evolve from this initial data 
lie in %. As discussed by Wald [2], since M is self-adjoint and positive, the 
answer is in the affirmative — further, the solutions evolving from initial 
data of compact support are defined everywhere, smooth, and bounded (in 
t) in the spacetime. Of course, we have obtained this result by choosing a 
particular boundary condition at the naked singularity. Horowitz and Marolf 
[I] have raised the issue of different choices of inequivalent boundary condi- 
tions leading to different evolution starting from a set of initial data. This is 
similar to quantum mechanical problems in which the Schroedinger operator 
is symmetric, but not self-adjoint. In such problems, there may be an in- 
finite number of self-adjoint extensions possible (i.e., inequivalent boundary 
conditions) each of which corresponds to a different spectrum. The choice 
of self-adjoint extension (boundary conditions) is dictated by the physics of 
the problem. In our case, we have done the same — only one unambigu- 
ous choice of boundary conditions implies regularity of the solution (so it 
remains a 'perturbation' and back-reaction effects can be ignored — other- 
wise the perturbative analysis in terms of the Klein-Gordon equation on the 
JNW background would break down). We have shown that for this choice 
of boundary condition, M is self-adjoint and positive. For other choices of 
self-adjoint boundary conditions involving the solution which is not regular 
at the naked singularity, the statements we have made on positivity of M, 
particularly after ( 1II.20|) are not valid. We do not know if those boundary 
conditions lead to positive self-adjoint extensions. Our choice of function 
space was the traditional space of square integrable functions with respect 
to a measure. For this very example, it has been mentioned in [5] that if the 
function space is an appropriate Sobolev space, M is essentially self-adjoint 
- and therefore there is a unique choice of boundary conditions at both 
endpoints. 

This completes the main arguments of this section, and our conclusion is 
that the JNW spacetime with a naked singularity is classically stable under 
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perturbations which are initially of compact support. 

Ill Discussion 

In this paper, we have shown that the JNW spacetime, which contains a 
naked singularity, is stable under scalar field perturbations which are regular 
at the singularity. This is a very surprising result, as we would have ex- 
pected naked singularities to be perturbatively unstable. Furthermore, the 
JNW spacetime is sourced by 'realistic' physical matter, a scalar field (unlike 
naked singularities in fluid models or dust, where the matter model breaks 
down at small scales). So the obvious generalization would be to consider 
gravitational perturbations of this spacetime and comment on its stability. 
If this spacetime turns out to be stable under gravitational perturbations as 
well, this would be indicative of the spacetime forming as a result of col- 
lapse for generic initial conditions. The question then would be if the cosmic 
censorship conjecture needs to be modified — it has been suggested, for ex- 
ample, that the physically relevant condition for a spacetime is not global 
hyperbolicity, but the well-posedness of all the probe field equations in the 
spacetime [28] . 

Another point of view could be that this naked singularity spacetime 
may have an instability in quantum gravity, even if it is stable classically. 
Thus, perhaps in those examples where naked singularities are classically 
stable, quantum gravity effects implement cosmic censorship. This surmise 
can be tested for this example — the question of stability of a spacetime in 
quantum gravity can be examined in the Euclidean path integral formulation 
[22] (also see [SU] and references therein). Thus, this stability result for the 
JNW spacetime opens up many tractable questions which will hopefully lead 
to a precise geometric/analytic criterion for stability of a naked singularity. 
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